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SU(6) proton wave function. 
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1. Introduction 

Quantum Chromodynamics predicts that single transverse spin asymmetries are suppressed 
in hard collisions, as a consequence of helicity conservation (chiral invariance) in the sub- 
process. These asymmetries indeed appear as interferences between helicity amplitudes 
which differ by one unit of helicity, therefore they vanish in the limit m quar k — > 0, or equiv- 
alently Q 2 — > oo (Q measures the hardness of the subprocess). Nevertheless, a number of 
high pt reactions persist in showing large asymmetries [1] . 

These facts do not invalidate QCD but mean that the approach to the asymptotic 
regime in p± is very slow, as regards polarization. However, in spite of their "nonasymp- 
totic" character, it is not unreasonable to think that the mechanisms of the asymmetries 
lie at the parton level. In other words, the asymmetries would be manifestations of quark 
transverse spin (or transversity) . Thus, we could extract information from them about 
the quark transversity distribution in the nucleon and/or the transversely polarized quark 
fragmentation. In this paper, we shall present a model for the spin asymmetry in the 
reaction 

p]+p^ir + X (1.1) 

which, unlike previous approaches [2,3], involves the transverse spin asymmetry of the 
polarized quark fragmentation [4,5], which hereafter will be refered to as the Collins effect. 

The paper is organized as follows: section 2 gives a very short review of the results of 
Fermilab E-704 experiment. In section 3, we explain how the Collins asymmetry can give 
rise to the observed single spin asymmetry in reaction (1.1) and deduce lower bounds on 
the transverse polarizations of the quarks in the proton, as well as on the size of the Collins 
effect. Section 4 presents a quantitative model based on string fragmentation and section 
5 gives the numerical results. Section 6 contains discussion of our results and conclusions. 



2. Main features of single spin asymmetry in inclusive pion production 

A strong polarization effect has been observed in the reaction (1.1) with 200 GeV trans- 
versely polarized projectile protons. The asymmetry is defined as 

A N (x F , Pl _)= a l~ a \ , (2.1) 

Assuming that j refers to the +y direction (vertical upwards), and the transverse momen- 
tum pj_ of the pion points towards the +x direction (pbeam is along the z axis). In other 
words positive An means that for upward polarization, the pions tend to go to the left. 
The most recent results, which we will consider in this paper, come from the Fermilab 
E-704 collaboration. They were published separately for two kinematical regions: 

- xf > or fragmentation region. The asymmetries have been measured for both 
charged and neutral pions [6-8]. 

- Central region, xf ~ [9], where the asymmetries were measured for neutral pions 
only. Large asymmetries for high p± in the central region had been observed previously 
also by other experimental groups [10] 

(xf is the Feynman variable 2p^ M / ^/s with p^ M being the longitudinal momentum of the 
pion in the CM frame). 
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In this paper we shall concentrate only on the first case. In this region the data 
show large asymmetry for all pions; positive for 7r + and n° and negative for tt~ . The 
asymmetries vary from about at xp ~ 0.2 to about +0.4, +0.15 and —0.4, for tt" 1 ", tv~ 
and ir~ respectively, at a;^ ~ 0.7 — 0.8. 



3. Hypothesis that E-704 asymmetry is due to Collins effect 

3. 1 Generalities from the parton model. 

In the "factorized" parton model (Fig. 1), the cross section for A+B — > n+X in the forward 
hemisphere is a kind of convolution of a parton distribution G q /A(x, <fj_), a parton scattering 

cross section & q +B^ q '+x = & q ^ q i and a parton fragmentation function D^j q i{z, h±) . In 
short-hand notations, 

(TA^tt ~ G q/A <g> Cfq^qi <g> D^/q, (3.1) 

[the factor o q ^ q i may be replaced by 5(q — q' ) (no scattering); this is the case of the Dual 
Parton Model]. Each factor in this equation may or may not depend on spin. Transverse 
polarization can act at three different levels: 

a) in a dependence of G q /a(x, q±_) on the azimuth of q± (Sivers effect [2]). 

b) in a single spin asymmetry in a q ^ q > (Szwed effect [3]). In this case, (but not necessarily 
in case a) , the quark q must inherit a part of the polarization of the proton. 

c) in a dependence of D n / q '(z,h±_) on the azimuth of h± (Collins effect [4,5]). Here, a 
transfer of polarization must occur not only from the proton to quark q but also from 
q to q' . 

The first mechanism implies nonzero internal angular momentum inside the proton, hence 
a strong departure from SU(6) model. As for mechanisms b) and c), SU(6) predicts : 



A N — 



--A" + 
2 



and 



A 11 
N 



< 



1 



(3.2) 



since most of the time, n + comes from a valence -u-quark, rr from a valence <i-quark of the 
projectile and the polarizations of u and d are respectively +| and — |. These predictions 

are in conflict with the E-704 results A 7 ^ ~ — A 1 ^ ~ 0.4 at large xp. To conclude, 
whichever mechanism we choose, SU(6) appears to be badly violated in G q / p (x, q±) . 



3.2 The Collins effect. 

According to Collins [4,5], the fragmentation function of a transversely polarized quark q 
takes the form 

D n/q (V q ,z,h ± ) =D n/q (z,h ± ){l + A 7r/q (z,h ± ) x \V q \ xsin[^)-^ ± )]} , (3.3) 



where V q is the quark polarization vector (\V q \ < 1), hj_ the pion transverse momentum 
relative to the quark momentum q and (p(a) the azimuth of any vector a about q. The 
factors after A can be replaced by \q x h±\~ 1 V q ■ (q x h±). Such a dependence is allowed 
by P- and T- invariance but has not yet been measured. 
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The Collins effect is the reciprocal of the Sivers effect. But Collins argued that the later 
is prohibited by time reversal invariance [4], while the former is not. As for mechanism b), 
it vanishes for massless quarks due to chiral symmetry: single spin asymmetry in q — > q' 
is not compatible with conservation of quark helicity. Therefore it should be small for 
hard or semi- hard scattering. In conclusion, among the sources of asymmetry a), b) and 
c) discussed above, we have a preference for the Collins effect illustrated by Fig. 1.** 

3.3 Consequence for the single spin asymmetry. 

Let us consider the hypothesis that E-704 asymmetry is due to the Collins effect. The 
polarized inclusive cross section reads 

Jp= ^2 J dxd2 ^ G q (x,q±) J dyd 2 k± G k (y,k±) x 

flavours of q, k, q' 

d(cos9)d(f - / dz d 2 h± D^/q^Vq', z,h±) S(p — zq' — h±) ; (3.4) 

dfl J 

the final quark transversity is given by 

V q , = K V heam A ^ x ^ d nn 0) . (3.5) 
G q {x,q±) 

1Z is the rotation about pbeam x q' which brings Pheam along q', 

A±G q (x,q±) =G^(x,q±)-Ggi(x,q±) (3.6) 

is the quark transversity distribution [12,13] in the proton polarized upwards, and Dnn(9) 
is the coefficient of spin transfer normal to the scattering plane in the subprocess. Formula 
(3.1) results from integration over the target parton variables y and k±. 

At large xf, the dominant quark flavours are q = q' = u for 7r + production, q = 
q' = d for 7r~ production. Furthermore, the hard scattering occurs predominantly at 
small 9 and Dnn{9) is close to unity, as in the case of t-channel one-gluon exchange 
[Dnn = —2su/(s 2 + u 2 )]. Assuming that A±G q /G q does not depend on <fj_, the results 
of E704 collaboration imply 

A±G u (x) x ^_ ^ Q 4 ^ (3 7) 



G u (x) 

A ± G d (x) 
G d (x) 



x A{z, h±) < about - 0.4 , (3.8) 



** We shall not discuss other approaches [11] not relying on the factorized parton de- 
scription (Eqs. 3.1 or 3.4). They are not necessarily in contradiction with the present 
one. 
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for xz ~ xp — 0.8. x means the most probable value of x. The inequalities take into 
account the fact that integration over q± and 9 allways dilutes the Collins effect. Thus, we 
get at least a lower bound of 0.4 separately for \A±G U /G U \ , \A±Gd/Gd\ at large x and 
\A(z, h±)\ at large z and for the most probable value of h±. 

4. Simple model of single spin asymmetry 

In order to make the conclusions of the previous section more quantitative, we per- 
formed a calculation in a simple model. 

We considered only the valence quarks of the projectile proton with distributions 
normalized as follows: 




[G q {x] = ^q va i(x) in conventional notations]. The quark is accompanied by a diquark 
carrying the fractional momentum 1 — x and whose distribution is 

Gt«(a,5l) = G d (l - x, -q ± ) ; Gyd(x,q±) = G u (l - x, -q±) . (4.2) 

We did not incorporate any hard or semihard scattering. Both 1Z and Dnn of Eq. (3.5) 
are equal 1. The two q — qq strings formed after the collision are parallel to the beam. We 
decay them recursively according to the simple Lund recipe [14]. We use the Standard 
Lund splitting function: 

f(z) = (l + C)(l-zf, (4.3) 

z being the fraction of the null plane momentum P + = P° + P 3 of the string carried by 
its leading hadron. f{z) = D rank=1 (z) corresponds to the production rate of the first-rank 
hadron (i.e. the one that contains the original quark spanning the string). This gives [14] 
for all ranks the inclusive density: 

D all ranks^ = ^ + ^ 1 ^ _ = \f( z ). (4.4) 

Z Z 

Thus, for all the other (subleading) hadrons originating from the string we get: 

D Tank ^ 2 (z) = (1 + C)i— ^(1 - z) c = f(z)^—^. (4.5) 

z z 

The transverse momenta of a quark and an antiquark of a pair created in the string balance 
each other (local compensation of the transverse momentum) and are distributed according 
to 

**l)JVx = ^«P (4.6) 

K \ K J 

k being the string tension. We took the intrinsic transverse momentum distribution 
in G q (x,q±) to be the same as for the tunnelling transverse momentum in the string: 
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G q (x,q±) = G q {x) p(q±). It yields an average intrinsic transverse momentum ^^ trinsic ^ = 
0.206 GeV. 

The polarization of the leading quark q is 



A ± G qo (a,g& ± ) 
G qo {x,q QA _) 



V qo {x)= Tr,T ; -V (4-7) 



and is taken to depend only on x. 

Each quark-antiquark pair created during string breaking is assumed to be in a 3 Pq 
state (vacuum quantum numbers) [15], i.e., with parallel polarizations. According to the 
Lund mechanism for inclusive A polarization [14], their polarizations are correlated to the 
transverse momentum of the antiquark q ± by 

*<-*' — ih-'-TT' (4 ' 8) 

where L is the classical orbital angular momentum of the qq pair and equals: 



K K 



(see Fig. 2 for a schematic explanation). 

In order that go and qi of Fig. 2 combine into a pion, they have to form a spin singlet 
state, the probability of which is 

I(l-7V7y, (4.10) 

in accordance with the projector on the singlet state \ — s(q ) ■ s(q~i)- The factor (4.10) 
causes the Collins effect: if qo in Fig. 2 is polarized upwards then qi and the pion which 
contains qi tends to go to the left-hand-side of the z direction. 

Vector mesons are ignored (accordingly, we omit the factor | in Eq. (4.10)). In our 
model the asymmetry for vector mesons would be of the opposite sign but three times 
smaller in magnitude when compared to that of the scalar ones. The observed asymmetry 
of the pions resulting from decays of the vector mesons would be even smaller, due to 
integration over decay ang lest. Anyway, at large xp, there are mostly direct pions. 

We do not introduce the Collins effect in subleading ranks or in the fragmentation of 
the diquark. For the fragmentation function (4.3), the probability that the detected pion 
is the leading one is equal to its momentum fraction z. Hence, the model can be a good 
approximation mostly for pions of high positive xf- We provide a short discussion of this 



t It has been shown however that the vector meson can also have a tensor polarization 
[16] which would result in the Collins effect for the decay products. We did not include 
this possibility. Another source of asymmetry [5] could be the interference between the 
vector meson and the nonresonating background. 
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point in the last section. We shall divide the sources of the pion production into three 
parts: 

1. the pion is the leading particle of the string spanned by the quark go- The Collins 
effect (and the asymmetry) appears only in this contribution. 

2. The pion is a sub leading particle of the string spanned by that quark, or 

3. it is a subleading particle of the string spanned by the diquark accompanying the 
quark qo- The leading particle of the string in this case is a baryon. Taking into 
account that xf ~ xz and that the transverse momentum of the pion is the sum of 
the transverse momenta of its constituents we get the production rates for all the three 
cases: 



da 



Otot 



dxpd 2 p_ 



rank=l 



= I dxdzd 2 q ± d 2 q ± G q (x,q ± ) a D rank=1 (z) 

quark q = u ,d^ 



x(l-P q . VSlS) P(q±) S(x F - xz) 5 2 (p ± - q ± - | ± ) (4.11) 



Ctot 



da 
dxpd 2 pj 



-i rank > 2 



quark 



E 

q=u,d ' 



dx dz G q (x) C2 D 



rank > 2 



(z) Pn(p±) S(X F -XZ) (4.12) 



Ctot 



da 



dxpd 2 p± 



rank > 2 « 

= dxdz c 3 D rank ^ 2 (z) P7T (p ± ) S(x F - xz) (4.13) 

diquark q=u ,d^ 



and the total production rate is the sum of the three: 



da 
dxpd 2 pj 



da 
dxpd 2 p± 



-i rank=l 



+ 



quark 



da 
dxpd 2 p± 



rank > 2 



+ 



quark 



da 
dxpd 2 p 



_L 



rank > 2 



diquark 



(4.14) 



Ptt(rl) is the convolution of p(q±.) and p(q±) and is also Gaussian but of twice larger 
variance^ . c\ , C2 and C3 are flavour factors and correspond to the probability that q and q 
match to form a pion of the appropriate charge. 

Since the production rate varies with the azimuthal angle of the pion momentum <p 
like in (3.3) then, in order to obtain the asymmetry at given values of p± and xp, we need 
to compare da(xF,P_i_, 4>) oruv at = and = tt: 



A n (x f ,p_l) = 



da(xp,P±, 0) — da(xF,P_i_, tt) 
da(x F ,P±, 0) + da(x F ,p±, n) 



(4.15) 



•t- We did not follow exactly the S(p — zq' — h±) prescription of (3.4) ; in other words 
we gave all the transverse momentum of the leading quark to the first rank particle. The 
resulting error is small at large xf- 
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which completes the calculation. 
5. Numerical results 

For the numerical calculations we parametrized the quark distributions as follows: 

G u (x) = 2G d (x) = ^x 1 / 2 (1-x). (5.1) 

We have used the string tension k = 0.17 GeV 2 and the parameter of the fragmentation 
function C = 0.3. In pair creation we have used the flavour abundances with the ratio 
u : d : s = 3 : 3 : 1, which determine the coefficients in Eqs (4.11-4.13) to be c\ = 3/7, 

C2 = cs = 9/49 for charged and C2 = cs = 18/49 for neutral pions§. 

In Fig. 3a we show our results compared to the data [7] of 0.7 < p± < 2.0 GeV. The 
dotted, dashed and full lines correspond to various dependences of the quark polarization 
V q on the momentum fraction x: constant, proportional to x and to x 2 respectively. 
Results obtained for all the three choices converge at xp = 1. The maximal values of the 
polarization, reached at a; = 1, are those motivated by the SU(6) wave function of the 
proton: V u = +2/3 and Vd = —1/3. 

One sees that the resulting asymmetry An strongly disagrees with the data for the 
negative pions. This confirms the conclusion of the Section 2. The measured asymmetries 
in tv~ production are, for xp >~ 0.5, equal in the absolute value but of the opposite sign 
to those of 7r + . This cannot be accounted for by SU(6) where the asymmetry of tt~ is 
negative but twice smaller than that of 7r + . 

The results compared to the data at lower transverse momenta, 0.2 < p± < 0.7 GeV, 
are shown in Fig. 3b. Here no disagreement with the data is seen. The data do not reach 
however as high values of xf as in the high p± interval. 

The behaviour of the measured asymmetries at high x f and p± encouraged us to try 
out a parametrization with maximal possible transverse polarizations of quarks one could 
choose in the proton: V u = 1 and Vd = — 1 at x = 1. The results are plotted in Figs. 
4a and 4b for both p± intervals. The constant polarizations V q (x) = const (dotted lines) 
result in too large asymmetry, at least at small and intermediate xp- This suggests that 
the transversely polarized quark densities fall down at small Bjorken x values. 

The full lines follow the data quite well. Only the measured asymmetries of n~ 
slightly differ from our results at small values of both xf and p±. One might conclude 
that A±G q (x)/G q (x) varies with x somwhere around x 2 but such conclusion is dangerous 
in scope of the simplicity of the model. 

In order to check how our results depend on the shape of the quark distribution G q (x) 
we did the calculation also for 

G^-x- 1 / 2 , (5.2) 

§ In this model, every uu or (id meson is considered as a 7r° (no rf) ; it gives a(n + ) + 
cr(n~) = o"(7r°), instead of 2 a(n°) as required by isospin. Nevertheless, Eq. (5.3) below 
remains true. 
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which differs strongly from (5.1) but gives better account of the leading baryonic effect 
(G^x) ~ (1 — x)~ x l 2 and the diquark tends to carry a substantial momentum fraction of 
the proton). The comparison is shown in Fig. 5. One sees that for the distribution (5.2) 
the asymmetry (dashed line) is slightly smaller at large xf but the difference is not large. 
The full line comes from Fig. 4 and corresponds to the distribution (5.1). 

In the parton model the 7r° asymmetry is just a combination of the tt + and tt~ ones : 

A i ox q(7r+) A N (n+)+a(n-) A N (ix~) 

A-NVK ) = 7 ~tt -, s • (5.o) 

Nevertheless, we show in Fig. 6 a comparison of our results to the E704 data [6] on n° 
production in pp and pp collisions. The curve obtained with V q {x) ~ x 2 (full line) agrees 
with the data also here. The lower p±_ bound in the 7r° case varies from 0.5 to 0.8 GeV 
depending on xp and was taken into account in our calculation. This is the reason of a 
bit wiggly shape of the lines. For the neutral pions the difference between SU(6) and the 
maximal polarizations of the quarks in the proton cancels and there is no difference in our 
predictions there. 

Finally, in Fig. 7 we plot the p± dependence of the asymmetry of tv°. The agreement 
with the data is also good. The shape of the curves reflects the Lund parametrization (4.8). 
Only the rise of the asymmetry at high p± and close to the central region of xp cannot 
be described by the model. We believe that this rise can result only from a combination 
of hard scattering with the Collins effect. The former has not yet been included in our 
calculation. 



6. Discussion and conclusions 

To summarize, we calculated the single transverse spin asymmetry in high-energy pp col- 
lisions in a simple model involving the Collins effect (asymmetry arising at the level of 
fragmentation of a quark into hadrons). We parametrized the Collins effect by the Lund 
mechanism of polarization in the coloured string model. 

We got good agreement with the data when we assummed that: 

a) The transverse polarization of the u and d quarks in the transversely polarized proton 
are close to unity but of the opposite sign (V u = P pro ton, 'Pd = —^proton) at momentum 
fraction x close to 1. 

b) The dependence of the quark transversity (or polarization) on the momentum fraction 
x is close to be proportional to x 2 . 

The conclusion a) is model-independent provided that the asymmetry arises in the 
quark fragmentation (Collins effect), which is a reasonable assumption as argued in sec- 
tion 2. 

The quark transversities we got : 

G u {x) G d (x) 



8 



are, in fact, not unreasonably large. Consider a covariant model of the baryon consisting 
of a quark and a bound spectator diquark [13,17] ; then 

°^ {X) = 16^ r ^ (J^) 2 ^ \u{q])Vu(p])\ 2 (6.2) 

00 1 ' diquark polarization 

where g(q 2 ) is the q — (g — B form factor, V = 1 for a scalar diquark, V = 757 • e for a 1 + 
diquark of polarization e M and 

Qm = xm B ~ T^ m <n ■ ( 6 - 3 ) 

J. X 

Formula (6.2) is similar to the covariant Weizsacker-Williams formula, but for a "spin ^ 
cloud"). Independently of g{q 2 ), this model predicts the following behaviours at x — > 1: 

- for a 1 + spectator diquark, helicity is fully transmitted [ ALG q (x)/G q (x) — > 1 ], 
transversity is fully reversed [ A±G q (x)/G q (x) — > —1 ]. In particular, Vd(x) — > — 1. 

- for a + diquark, Aj_(jr g (a;) and G q +(x) coincide and, for g(q 2 ) decreasing faster than 
q~ 2 , they exceed ^G q (x) as x — > 1 

Thus, a dominance of the scalar spectator for the it quark and the pseudo-vector one for 
the d at x ~ 1 could lead to the large opposite transversities as in (6.1). 

The conclusion b), related to the xp dependence, is model-dependent and cannot 
be taken too seriously. One needs a good parametrization of the Collins effect before 
one can deduce the x dependence of the quark transversity. Our parametrization is an 
approximation which should work only at reasonably high values of both xp and pt- We 
took into account the Collins effect only for the first-rank (leading) hadrons, wherefrom 
A oc z in (3.3). The second-rank hadrons have the asymmetry of the opposite sign as 
compared to the first-rank ones. More generally, the subsequent ranks are asymmetric 
in the opposite way to each other (as required also by local compensation of transverse 
momentum). This should cause a faster decrease of A at lower z values, where the higher- 
rank hadrons are more important. Unfortunatly this feature was not possible to include in 
our simple semi-analytical calculation, since the yields of rank-2 (and higher) hadrons do 
not have simple analytical forms. The contribution of vector mesons also should reduce A 
at lower z. Assuming x ~ z ~ yfxp, a steeper A (for instance oc z 2 ) would imply a flatter 
V q {x) (for instance oc x). 

In our model the high p± hadrons originate from the tail of the Gaussian distribution of 
the intrinsic and tunneling transverse momenta. Incorporating the hard scattering should 
not change the results very much at large p±. There will be some smearing of h±, hence a 
reduction of the asymmetry, but not too severe because of the "trigger bias" effect: in most 
high p± events, the contributions of string decay and of hard scattering to the transverse 
momentum are rather large and point in approximately the same direction, as in Fig. 1 
(in order to sum up to the large p± of the pion). At small p±, the trigger bias effect would 
work less and this could improve the agreement with the data in that region (Fig. 4b). 

One main conclusion of this paper is that single spin asymmetry may be the first 
experimental indication for the existence of the Collins effect. A more detailed experiment 
would be usefull to select between this and alternative explanations. Besides its theoretical 
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interest, the Collins effect may be the most efficient "quark polarimeter" necessary for the 
measurements of the transversity distributions in the nucleons [5,18]. We hope that this 
effect will soon be tested directly, for instance in the azimutal correlation of two pion pairs 
from opposite quark jets in e + e~ annihilation. 
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Figure captions 

Fig. 1 Inclusive pion production. Two events (a) and (b), symmetric with respect to the yz 
plane, are represented. Without polarization, they would have the same probability. 
In the polarized case, the Collins effect favours the case (a). 

Fig. 2 Production of the leading pion in a string drawn by a transversely spinning quark. 

Fig. 3 Single spin asymmetry measured by E704 collaboration. The curves are our model 
results calculated with quark transverse polarizations at x = 1 as in SU(6) wave 
function of the proton (V u = 2/3 and Vd = —1/3). The data are from Ref. [7]. 

Fig. 4 The asymmetry calculated under an assumption that the quark polarizations at x = 1 
are: V u = +1 and Vd = —1- The data are as in Fig. 3. 

Fig. 5 Comparison of two various quark distributions. The full lines and the data are as in 
Fig. 4. 

Fig. 6 The asymmetry of 7r°'s. Here there is no difference between the predictions of SU(6) 
and the maximal polarizations. The quark distributions are as in Figs 3 and 4. Data 
come from Ref. [6]. 

Fig. 7 p± dependence of the ir° asymmetry in two intervals of the Feynman variable xp. 
Data are from Ref. [8] . 
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